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however, not been studied at all. In this paper we look at three classes 
of stable random fields, and obtain asymptotic formulae for the mean 
values of various geometric characteristics of their excursion sets over 
high levels. 

While the formulae are asymptotic, they contain enough information 
to show that not only do stable random fields exhibit geometric be- 
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1. Introduction 

We are interested in the structure of the sample paths of certain smooth 
stable random fields / : R w — > R, N > 1. We shall study these through 
certain geometric properties of their excursion sets 

A u = A u (f,M) = {i E M : f(t) > u} , (1) 

where M C M N and u G R. 

Excursion sets have been widely studied for Gaussian and Gaussian re- 
lated random fields. Their applications appear in disciplines as widespread 
as astrophysics and medical imaging, where they have been also been used 
in a variety of hypothesis testing situations A good introductory reference 
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to the applications is still Keith Worsley's exposition [9] although |4j, when 
ready, will have a lot more detail. On the more theoretical side, where ex- 
cursion sets are seen to generate an elegant geometric structure, our basic 
reference will be the recent monograph [3 . 

The Gaussian and Gaussian related scenarios allow for the development of 
explicit formulae for the expectations of many of the geometrical quantifiers 
of excursion sets. Unfortunately, one cannot expect that the same will occur 
in the stable case, for which our basic reference will be [8]. Here, explicit 
formulae for even the marginal densities of / are unknown, although much 
is known about their asymptotics. Indeed, when we began this research we 
expected to be able to find little beyond some asymptotic formulae relating 
to the excursion sets generated by stable fields that might mimic the Gaus- 
sian ones, much as we did for level crossings of stable processes on the real 
line, in [I] and [2]. 

What we found turned out to be far more interesting. It is well known 
that the structure of stable processes is far more complicated than that 
of their Gaussian counterparts. For example, whereas in the Gaussian case 
many stationary processes have both a moving average representation (with 
respect to white noise) as well as a harmonisable representation, in the sta- 
ble case moving average and harmonisable processes belong to quite distinct 
families. These differences are well understood in terms of mathematical 
structure, but it turns out that the sample path distinctions between differ- 
ent classes of stable processes become significantly highlighted by looking 
at the excursion sets that they generate, primarily in the multi-dimensional 
setting. 

Thus, this paper has two aims. One is to provide explicit, albeit asymp- 
totic, formulae for the expectations of the Euler characteristics (defined be- 
low) and other geometric quantifiers of excursion sets generated by stable 
random fields. These, we believe, will find immediate application in random 
field modelling using these processes. The second is to better understand 
the differences between various stable random fields via their excursion sets. 

The remainder of the paper is structured as follows. In the next sec- 
tion we shall define the Euler characteristic and related Lipschitz-Killing 
curvatures. In Section [3] we give a description of the main result from the 
Gaussian theory, which is a precise formula for the expected value of the 
Lipschitz-Killing curvatures of Gaussian excursion sets. With this in hand, in 
Sections |4]^6] we start with the new results, for sub-Gaussian, harmonisable, 
and concatenated-harmonisable random fields, deriving asymptotic formulas 
for the expected values of the Euler characteristics of their excursion sets. 
In Section [7] we show how to lift these results to all the Lipschitz-Killing 
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curvatures, at least for isotropic fields. A technical appendix completes the 
paper. 

2. Euler characteristics and Lipschitz-Killing curvatures 

Throughout this paper our parameter sets M will be taken to be convex 
subets ofR N , and, for much of it, we shall restrict ourselves to A^-dimensional 
rectangles. In these cases we shall always write T rather than M, where T 
is given by 



For convex sets in M. , and for excursion sets of smooth functions defined 
over them, there are N + l functionals which describe their geometry. These 
are known under a variety of names, including Minkowski functionals, in- 
trinsic volumes, quermassintegrals, and Lipschitz-Killing curvatures, being 
related to one another by differences in the way they are ordered or normal- 
ized. For consistency with [3], which we shall use heavily when we wish to 
cite a result without proof, we shall work with Lipschitz-Killing curvatures. 
Perhaps the easiest way to define these is via Steiner's formula, a classic 
result of integral geometry. 

To state Steiner's formula, let M be a convex set of dimension N, sitting 
in ~K N , where N' > N. (e.g. M is a one-dimensional curve in M 2 .) The tube 
of radius p around M is defined to be 



With Atv denoting Lebesgue measure in K , Steiner's formula states that 
there is an exact polynomial expansion of order iV for the \n' measure of 
Tube(M, p) given by 



N 



(2) 



3=1 



Tube(M, p) = {x e R N ' : inf \x-y\< p}. 



N 



\ N , (Tube(M,p)) = Y / u N >-jP N '~ j C j (M) 

j=0 



(3) 



where 



•i/2 
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is the volume of the unit ball in W . The numbers Cj{M) are the Lipschitz- 
Killing curvatures of M. 

The Cj scale nicely, in the sense that, for A > 0, Cj(XM) = A J £j(M). As 
is obvious from Cn{M) measures the volume of M, Cn-\(M) is related 
to its surface measure, etc. The last one, £q(M), is the Euler characteristic 
of M which, since we shall use it very often, we also denote by <p(M). Of all 
of the Cj(M), the Euler characteristic is the only one that does not change 
under smooth deformations of M. If M is one dimensional, then <p(M) counts 
the number of connected components in M. If M is two dimensional, then 
it counts the number of connected components minus the number of holes. 
In three dimensions, it counts the number of connected components, minus 
the number of 'handles', plus the number of holes. 

We note, for later usage, that the Lipschitz-Killing curvatures of the Af- 
reet angle ([2]) are given by 

N 

C J (T) = C j (l[[0,T i })=Y / T il ...T l] , (4) 

i=l 

where the sum is taken over the ( . ) distinct choices of subscripts i%, . . . , ij. 

There is another way to write Q, which will be useful later on. Let 
Oj = Oj (T) denote the collection of the ( ■ ) j-dimensional facets of T which 
contain the origin. (Thus, for example, Ojy is T itself, while 0\ contains the 
N one-dimensional edges of T lying on the positive axes e\, . . . , ejy of 1^.) 
Furthermore, if J is a facet in Oj, let \J\ denote its j-dimensional Lebesgue 
measure. Then it is immediate from Q that 

C j (T)=J2\J\- (5) 

The Lipschitz-Killing curvatures play a central role in much of integral 
geometry, but for the moment we shall note only one of their properties, 
known as Hadwiger's theorem [6]. Suppose that we have a functional ip on 
compact convex sets that is additive, in the sense that, if A, B and A\J B 
are compact convex, then 

tp(A U B) = ip(A) + tp(B) - ip(A n B). (6) 

If it is also true that ip is invariant under rigid motions and continuous in 
the Hausdorff metric, then there are (^-dependent) constants such that 

N 

^{A) = Y J c j C j {A). (7) 

3=0 
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Thus, studying intrinsic volumes is equivalent to studying a far wider class 
of functional on sets. Our aim is to study the Lipschitz-Killing curvatures 
of excursion sets. 

3. Gaussian excursion set geometry 

In this section we want to summarise some results about the excursion sets 
of Gaussian random fields. There are two reasons for bringing these. The 
first is that it gives us a basis to which to compare the results of this paper 
for stable random fields, and the second is that, in all the cases that we 
shall consider in this paper, the proof for the stable case follows from the 
Gaussian one and a conditioning argument. 

To state the main result for the Gaussian case, we need a little notation, 
for which we now assume that / is a mean zero, stationary, Gaussian random 
field on R N with constant variance a 2 . Assuming that / is also almost surely 
C 2 , we define the second order spectral moments 



With {ei, . . . , e^} denoting the positive axes of R , suppose that J € Oj is 
a j-dimensional facet of a rectangle T. We write Aj for the matrix 



Note that if / is isotropic then there is a constant, which we shall write as 
A2, such that 




Aj = {Xij : i, j € cr(J)} , 



(8) 



where 



a(J)±{j: ej n(J\{O})^0}. 




(9) 



Next, we need the Hermite polynomials 



[n/2\ 

H n (x) = n\ ]T 



(-iy x n - 2j 

j\(n-2j)\2i 



n > 0, iGR, 



3=0 
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where [a\ is the largest integer less than or equal to a, and, for notational 
convenience, we define 

where 




We also adopt the notation 

p n (u) = (2vr)-(" +1 )/ 2 e"" 2 / 2 , n > 0. 

The following is a combination of Theorem 11.7.2 and the discussion in 
Section 11.8 of [3]. 

Theorem 3.1. Let f be a zero mean, stationary Gaussian field on a N- 
rectangle T with variance a 2 and a.s. C 2 sample paths, and such that the 
joint distribution of f and its first and second derivatives at each point t £ T 
is non-degenerate. Suppose that the joint modulus of continuity uj(rj) of all 
the second order partial derivatives of f satisfies 

P {u{rj) > e] = o 0q N ) as T] | (10) 

for all e > 0. Then the mean value of the Euler characteristic of its excursion 
set is given by 

n =ojee>„ ° vcr/ 

Furthermore, if f is isotropic and M compact and convex, then, for all 
0<J<N, 

N-j 

E{C j (A u (f,M))}=Y, 

n=0 

where A2 is as in Q and 

~N 

j 



C n+j (M)p n {^jl^J , (12) 



N 



J I LO N -j Uj 
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In fact, results of this nature hold in much wider generality, when M is a 
general stratified manifold M and / is neither stationary not isotropic. We 
refer the interested reader to Chapters 12 and 13 of [3]. 

One observation that follows from (11) and (12) comes by rewriting the 
sums as power series in u, from which one immediately sees that the leading 
order terms, of orders u N ~ 1 and u N ~^~ l respectively, are associated with 
the volume of T or M, while lower order terms are associated the other 
Lipschitz-Killing curvatures. This observation will be important when it 
comes to understanding the stable case. 

Finally, although it will not be important in what follows, we note that 
the condition (10) is an extremely mild one for C 2 Gaussian processes, and 
is easily checked from the covariance function of the process. 



4. Sub-Gaussian fields 

With the Gaussian case behind us, we shall now look at what is probably 
the simplest of all stable random fields, the sub-Gaussian ones. Despite their 
simplicity, we shall see that their behavior is already very different from the 
Gaussian case. 

To define these processes, we let g be a Gaussian random field on M, and, 
for some a S (0,2), let X be a ^^(ctq,, 1, 0) random variable independent 
of g (see [H] for notation) , where 

a a = cos(ira/4) 2 / a . 

Thus X is a positive strictly a/2-stable random variable with Laplace trans- 
form 

E{e- tx } = e- ta/2 , t>0. 

Taking X independent of g and setting 

f(t) = X 1 / 2 g(t). 

defines a sub-Gaussian random field. 

To state our first result we need some notation. For any functions a, b : 
R^Kwe write 



a x b 
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Theorem 4.1. Let g be a zero mean, stationary Gaussian field on the N- 
rectangle T, satisfying the conditions of Theorem \3.1\ Denote its variance 
by (Tg and the matrices of its second order spectral moments by Aj, as in 
Q. Let f, as above, be the sub-Gaussian field X l l 2 g. Then 

E{^(i»(/,r))}x M - fl k + E^n £ \J\\Aj\ l ' 2 Y ( i3) 

V n=l JeCV / 
where we write \Aj\ for |det(Aj)|, 

o-l+a/2 n .a r ( a+l \ 

_ a g L { 2 ) 

Ao " v^T(l " f ) ' 

and, for n > 1, 

a2 -l+«/2 CT a r («±l) (n _ x) , L(n-1)/2J (_lpT (°±"=l=3i) 

K ™ = ~ r(l - f) vr (n+i)/2 CJ n £ 2 2 i+ 1 j!(n-l-2j)!' 



If g ( and so f) is isotropic, with second spectral moment A2 (cf. and M 
is a compact convex domain, then 

N 

E {99 (A u (f, M))} x u- a ]T K n \ n 2 /2 C n {M). (14) 

n=0 

Under isotropy, a corresponding result holds for the expected Lipschitz- 
Killing curvatures of other orders as well. We shall look at this later, in 
Section 

Before proving this theorem, we shall take a moment to see how very dif- 
ferent it is from the purely Gaussian case, despite the fact that the excursion 
sets of / and g are very simply related by the fact that 

Mf,T) = A uX . 1/2 (g,T). 

Consider the case N = 1, when the theorem relates to a / defined over the 
interval [0, T]. Then 

ip{Mf, [o,r])) = i /(0) > u + c u (/,T) J 

where C u (f, T) is the number of upcrossings of the level u by / in [0, T\. Tak- 
ing expectations, we see that the two terms that appear in ( 13 ) correspond 
to (asymptotics for) E{lj( )>«} = P{/(0) — u } an d 



E{C n (/,T)}xu 



a 2 -iW2r(i + |)(A 11 )V2 r 

(1 _ |) a o-*w ' 
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where the asymptotics here come from either substitution in ( 13 ) or Theorem 
3.2 of [2j, which studies one dimensional level crossings. Note that both terms 
- i.e. the "boundary" and "interior" terms - have the same asymptotics, of 
the form u~ a . This is also true when looking at the case of general N in 



( 13 ), in that facets of T of all dimensions contribute to the asymptotics, and 
this is probably the most interesting aspect of the result. 

Recall that in the Gaussian case we saw that in an expansion of the mean 
Euler characteristic the leading term involved only the volume of T, with 
the surface area affecting only the second and later terms, etc. In the sub- 



Gaussian case, however, it is clear from (13 1 that the full geometry of T 
affects the first term of any such expansion. 

A heuristic explanation for this is easy to find. In the Gaussian case, 
if the level u is high, the excursion set will, with high probability, con- 
tain only a small set which will be unlikely to intersect the boundary of T. 
Hence, only volume terms appear in the highest order term when expanding 
E {</? (A u (f, T))}. However, in the notation of the theorem, an excursion set 
of the sub- Gaussian / at the level u has the same geometry as an excursion 
set of the Gaussian g at the level u/\/X. Although u may be large, the most 
likely reason for / to reach this level is that X also be large and that uj\[X 
is roughly O(l). This being the case, A u ^(g,T) is a "typical" rather than 
"rare" excursion set for g, and so has a reasonable probability of meeting the 
lower dimensional facets of T. Thus it is not surprising that they contribute 



to ([13]). 

Despite the (hopefully) convincing tone of these heuristics, the proof fol- 
lows a different, and purely analytic, route. 



Proof of Theorem 4.1 . We shall start with the general, non-isotropic case, 
and T a rectangle. 

As noted above, A u (f,T) = A uX -i/ 2 (g,T). Thus it is immediate that 
(p(A u (f, T)) is well defined since this the same is true of ip(A u (g, T)), for 
every u. Conditioning on X, we would now like to claim that 

E{<p(A u (f,T))} = E{E{<p(A u (f,T))}\X}, (15) 



and then use the Gaussian Theorem |3.1| to compute the inner expectation. 

However, to justify this we need to establish two facts. The first is that 
the conditioned process, f\X, which is clearly Gaussian, satisfies all the 
conditions of Theorem |3.1| This is trivial, since g is assumed to satisfy these 
conditions and f\X is no more than a constant multiple of g. 



The trickier problem is that to apply the iterated expectation in (15) we 



need to know, a priori, that the absolute moment E {\ip (A u (f, T))\} is finite. 
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Since this is technical we shall leave it to Lemma 8.1 in the appendix, and 
for the moment progress assuming that it is true. 



Then, applying Theorem 3.1 to (15), we have 
E{<p(A u (f,T))} 

= E{E{<p(A u/V x(g,T)) \ X}} 



N 



U 



n=0 



v \J\ \Aj\^ 



To evaluate this triple sum (the third sum appears implicitly in the Hermite 
polynomials) we need to consider typical terms of the form 

E{u k {a 2 g X)- k ' 2 e- u2 l 2 ^ x } = (^j E {x^/V^X j 
and the one atypical term, coming from n = and H—i, of the form 




The asymptotics of these expressions is covered in Lemma 4.2 which will be 
crucial to most of the computations of this section. Once we prove Lemma 



4.2 (13) is the consequence of a little algebra, and so the first part of the 



theorem is established. 

As for the isotropic case, note that taking M = T, (14) follows immedi- 



ately from ( |13[ ) on noting (|5) and the fact that detAj = Af im J . To establish 
the argument for general e convex compact M we use this fact and Had- 
wiger 's result @. 

To this end, define a functional ip on convex compact sets M C M. N by 
setting 

iP(M)^E{ip(A u (f,M))}. 

It is immediate that ip is additive (in the sense of ([6])) and, using the same 
general arguments as in the appendix, by bounding ip(A u (f,M)) by the 
number of critical points of / over M, that ip is also continuous over convex 
compact sets. Furthermore, if g is isotropic, it follows that ifi is also invariant 
under rigid motions. Thus Hadwiger's result applies and we have that there 
exist constants Cj such that, for convex M, 

N 

E{<p(A u (f,M))} = J2cj£j(M). 

j=0 
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Since the Cj are not dependent on M, the fact that ( 14 ) holds for rectangles 
defines them, and this fact and the above equation implies that (14) holds 
for convex M as well, as required. □ 



Lemma 4.2. Let X be as in Theorem 4-1 Then, for any (3 > —a, 



lim/ + %(l^/ 2 e - 2 ^^) (i 6 ) 



a + (3 



Furthermore, 



Proof. The limit (16) is a trivial consequence of Lemma 2.2 of pp. As for 



(17), let G be a standard normal variable. Then 



iE{|G|"}r{x> u 2 /^} 



the last line following from a classic result of [5]. 

Since E{|G| Q } = 2 a l 2 ^- 1 ' 2 T{^) and F{X > v} x C a/2 a% /2 v- a / 2 (see 
(3.7.2) and (1.2.8) of [8 d , respectively), the result now follows. □ 



5. Harmonisable fields 

The sub-Gaussian random fields of the previous section provide an inter- 
esting class of processes, in that they show that even a relatively minor 
perturbation of the Gaussian scenario leads to quite different behavior of 
the excursion sets. However, they do not represent a particularly rich class 
of stable fields. A much richer class of stable fields is given by the station- 
ary, symmetric, a-stable (SaS) harmonisable ones. These are random fields 
possessing a spectral type representation of the form 



/(i)=Re{ f e^Zidu)}, 



(18) 
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where Z is a complex, SaS, Borel random measure on M. N with finite control 
measure //. (See [8] for details of this and the following representation.) 



While (18) may explain from where the terminology comes, there is an 



alternative representation that will be much more useful for us, and which 
is given by 

oo 

f{t) = (C a b- Vo) 1/Q E T ~k l ' a ( G k ] cos (^) + sm(^ fc )) , (19) 
fc=l 

H N and the product tu^ is actually the inner product Yli=i ^W^fcW- 
, i = 1, 2, are independent sequences of i.i.d. standard normal vari- 
ables, {r^} is the sequence of arrival times of a unit rate Poisson process, 
{u>k} is a sequence of i.i.d. 1^ valued random variables with probability 
measure fj,(-)//j,o where 



where i € 1 
The{Gf} 



/x = Af(K 



(20) 



and /i is the control measure for (18). The four sequences are independent 



of one another. The constants C a and b a are given by 



and 



: sin x dx 



-i 



'(r(l - a)cos(vra/2))- 1 if a / 1, 
2/tt ifa = l, 



6 a = 2°/ 2 r(i+ f). 



An important consequence of the representation (19) is that if we condi- 



tion on the sequences {r^} and {uik} then the conditioned field is stationary 
Gaussian. This will enable us, as in the sub-Gaussian case, to use conditional 
Gaussian arguments to prove the following result, in which we establish the 
asymptotics of the expected Euler characteristic of the excursion sets of the 
real harmonisable stable fields. For simplicity we restrict ourselves to the 
case of a compactly supported control measure fi with a bounded density, 
but we expect the result remain true in greater generality. 



Theorem 5.1. Let f be a harmonisable, SaS, random field as in (18) or 



(19), defined on the N -rectangle To/Q. Assume that the control measure \x 
has compact support and a bounded density with respect to Lebesgue measure. 
Then 



E{<p(A u (f,T))}y< U - a C a no 



' 2 l-a/2 r (ifca) 1 



N 



N TA, (21) 
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where fio was defined at (20) and the ftj, j = 1, . . . , N are the (normalised) 
moments 

M.MN^/N^- (22) 
J no 

If, furthermore, ft is rotationally invariant (so that f is isotropic) and M is 
a compact convex domain, then 

(23) 

2l-a/2p ( l+a N 



E{cp(A u (f,T))} x u- a C af i \ ^ 2 J C (M) + ^Ci(M) | , 



where fil *s an U of the (equivalent) moments given by (22). 



Note, once again, how different this result is to the corresponding Gaus- 



sian one, and even to the sub-Gaussian one. Comparing, for example, (23) 



with its Gaussian counterpart ( 12 ) (take j = there) we see that while the 
leading term in the Gaussian case comes from the volume term Cn(M), in 
the harmonisable stable case it comes from the two lowest Lipschitz-Killing 
curvatures, Co(M) and Ci(M). 

To see why this should be the case, we shall postpone the rather technical 
proof of the theorem for a moment, take a moment to describe the principles 
involved, and then use them to obtain a heuristic proof of the theorem. 



Consider the representation (19) for harmonisable stable fields. From this, 



one can argue that, conditional on / reaching a high level u, the first term, 
with the coefficient T 1 1 , will dominate not only all the other summands, 
but in fact their sum. (Formulating this properly, and then establishing it, 
is basically the main part of the technical proof.) 

This being the case, / will tend, at high levels, to look like a cosine 
function, with random height, frequency, and direction. Thus, for example, 
in M 2 , a typical high level excursion set will look like that in the square of 
Figure [TJ a sequence of strips looking like the tops of cosine waves, with 
small perturbations due to the terms in the representation other than the 
dominant one. 

To count how many such strips there are, one needs only to look at the 
boundary of the square. In fact, if we increase the size of the square, it 
is clear that the number of strips (and so the Euler characteristic) grows 
proportionately to the length of the edges, and not to the area. Thus it 



should no longer be surprising that only £q(M) and £i(M) appear in (23). 
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(The Cq{M) term arises to 'catch' regions such as those in the lower left 
corner of Figure [TJ) 

There is also another argument that will give us not only the results of 



Theorem 5.1 but even a little more. Again, it is only heuristic, but it is both 
elegant and simple, and since it may also be applicable to other problems it 
is worth the space we shall devote to it. Its approach is via integral geometry. 
We start with the assumption that it is only the first term of the sum 



(19) that will be important, and so (ignoring multiplicative constants) look 
at the random field 

f{t) = r^ 1/a (GS 1} cos(^i) + GS 2) sin(twi)) 

over a compact, convex M. 

For a given direction to E R , let denote projection onto the line 
containing the origin and u and let 

M u = {P w t : t£M} 

be the projection of M onto this line. Conditioning on lu, the Euler character- 
istic of the excursion set A u (f, M) is the same as the number of upcrossings 
of the level u by the one-dimensional cosine wave on , plus one if / > u 
at the boundary point of M u closest to the origin. 

If we now average over the Gaussian variables Gi and then it is easy 
to check that, conditioned on uj, the expected Euler characteristic, should 
be proportional to 



F{/(0)>u} + u- a |M||Tj xn- Q (A;i + k 2 \\uj\\\T L 



for some constants k\ and k 2 which we shall not worry about. Averaging 
over u we find that 



/ ( kl + k 2 \\u\\\T^\) ^> = h + k 2 f \T Q 
Jr n Mo Jr n Mo 
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Now assume that the measure \x is rotationally symmetric. Then the so- 
called projection theorem of integral geometry (cf. [7j, §7.4) gives us that 



the integral here is proportional to C\(M). This gives (23), once we compute 
the constants. 



Indeed, we can go further than (23), dropping the isotropy assumption. 
In general, the quantity ||u;|||M w | can be expressed as 



where 



sup(o;,i) — inf {u,t) = /im(^) — %/(— w) 



/im(^) = sup(u, t) 



(24) 



(25) 



is the support function [7] of the convex body M. (Usually, the support 
function is considered as a function on the unit sphere S(M N ), but (25) is 
well defined as a function on M. N .) 

Following through with the constants, this argument would give that, for 
compact convex bodies M, 



E{<p(A u (f,M))} 
vT a CaHa 



2 l-«/2 r (!±«) 



+ 



1 

2^ 



(h M (u) - h M (-uj)) 



Mo 



(26) 

To see how this works for rectangles T, note that in this case the difference 



of the support functions appearing in (24) is just J2jLi \ u j\Tj- Substituting 
this into (26) gives us back (21). 



We shall return to integral geometric arguments later, in Section [7j but 
for the moment we leave heuristics and geometric arguments and give the 
promised technical proof. 



Proof of Theorem |5.1[ We shall only give a proof of ( 2 1 ) , the result for the 
rectangular parameter space T. The extension to compact, convex parameter 
spaces follows from Hadwiger's representation of additive functionals, as in 



the proof of Theorem 4.1 



The proof will follow the general lines of that of Theorem 4.1 in that it 
begins with a conditioning argument from which asymptotics are computed. 
The harmonisable case, however, is somewhat more complicated. 



We begin with the representation (19). If we condition on the sequences 



{T/j} and {uJk}, it is immediate that the conditioned field is a stationary 
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Gaussian field on M with mean zero and covariance function 



-2/a 

J H l ; - la Z^'-k 
k=l 

where 



^) = 7«E r fc 2/Qcos (^)> 

k=l 



7« = («Vo) i/a (27) 

Therefore it has variance 

oo 

5 2 =7^E^ /a (28) 

k=l 

and second spectral moments 



X ij = 7a $3 r fc w *(O w fc(i)- 
fc=i 



Our next step will be to apply the Gaussian Theorem 3.1 to the condi- 
tioned harmonisable process, and then use the fact that 

E {<p(A u (f,T)} = E{E {if (A u (f,T))\T k , co k , k > l}} . 

As in the sub-Gaussian case, the iterated expectation requires justification, 



which will be provided only later in Lemma 8.4 

In order the apply the Gaussian result, we need first to verify three condi- 
tions: that the conditioned process is a.s. C 2 , that the joint distributions of 
the various conditioned derivatives are non-degenerate, and that the condi- 



tion (4.1 ) on the moduli of continuity of the conditioned process is satisfied. 

We tackle the first of these first, by showing that /, itself, is a.s. C 2 . From 
this fact and Fubini's theorem, the same will be true of the conditioned 
processes. 

Note that the assumption of compact support for the control measure // 
implies, by Corollary 11.7.5 of [3], that / is absolutely continuous. Write 
fi and fij for the various first and second order partial derivatives of /. 



Exploiting the representation (18), it is easy to check that as versions of its 



partial derivatives we can take the random fields 

fj(t)= Re{ / ^iu(j)Z(du;)}. 



Since these are, again, harmonisable fields with compactly supported control 
measures they are all continuous. From this it follows that / is C . Applying 
the same argument to the derivatives of the /,■ gives that / is C 2 . 
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We now turn to the issue of the non-degeneracy of joint distribution of 
the various derivates. Applying the representation (19), we have 

oo 

h(t) = (c Q 6-Vo) 1/a £ r* 1/a W(t) Mt"k) + 4 2) cos(^ fc )) , 



k=i 



and 



(29) 



kit) 



oo 

-{C a b~ x no) Xfa r fc 1/a w fc («Vfc(i) (G^ sm(tu k ) + G ( k z> cos(tw fc )j , 

k=l 



i,j = 1,...,N. In particular, the joint distribution of / and its partial 
derivatives at time zero (up to a positive multiplicative constant) is that of 
the random vector 



(i) 



.fe=i 



fe=i 



/c=i 



Since we have assumed that the distribution of the random vectors ujk has 
a no n- vanishing absolutely continuous component, this is, with probability 

I, a non-degenerate Gaussian vector. 

It remains to check that condition (10) on the moduli of continuity of the 
second order derivatives of the conditioned process is satisfied. By Corollary 

II. 3.2 of [3], writing Cf i {t) for the covariance function of the second order 
partial derivatives /y, this condition will be satisfied if, for small enough t, 
and some finite K, rj > 0, 



max 

hi 



< K\]n\t\ 



-(!+??) 



(30) 



However, using the representation (29) for the it is easy to compute an 
expression for the CV, and, using the compact support of the uJk(i), to see 
from this that (30) is not only satisfied, but that the bound on the right 
hand side can be taken of order Itl 2 . 



With all the conditions checked, we can now apply Theorem 3.1 to the 
conditioned harmonisable process, from which it follows that 



E{<p(A u (f,T))} 



(31) 



N 



E I r- J 2 - V V |J||Aj|1/2 



H, 



n-l 



II 



a 



II 
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where the expectation in the second line is over the and Uk, and for 
J £ O n the nxn matrix Aj bears the same relation to the \ij that Aj does 



to the Xij. 



The term E (u/a)} can be handled much as in the proof of Theorem 



4.1 using (17) to show that 



2 1 -°/ 2 c Q/ x r (*±g) 

y/7rb a 



(32) 



As for the other terms, it is clear, expanding the Hermite polynomials, that 
we need to study the asymptotics of terms of the form 

{|A |l/2„,n-l-27 „ 1 

for < j < L^J • 

To do this, we need a little notation. For k > 1 and a facet J of dimension 
n define the nxn matrices Wj. = Wfc(J) by setting 



Wk(i,j) = w fc (i)w fc (j), i,j £ a{J). 



Then (32) can be rewritten as 



(33) 
(34) 



7, 



1/2+% , 



1/2 



U 



n-l-2j 



-u 



^oo r -2/o\(2«-l-2j)/2 



exp 



Our first step in handling this expectation will be to truncate the sum in 
the numerator. Note that, for any two nxn matrices A and B, the standard 
expansion of a determinant shows that 

n 

det(4 + B) < det(A) + £ A][ r * _w, '[| J B]| m , 

m=l 

for some combinatorial constants C nm that we allow to change from line to 
line and \\A\\ |. Furthermore, since y/x + y < \fx + -y/y, we can 

apply this to the numerator of (34) to see that there is a constant C such 
that 



(det(Er fc - 2/ ^ fc )) 1/2 -(det(r- / ^ 1 



1/2 



(35) 



m=l 



v-2/eti 



(n-m)/2 



E^ 2/Q ll^l 



fc=2 



m/2 
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We now claim that, in view of the above inequality, the expectation in (34) 
differs from 



det (i 



-2/a 



1/2 



,n-l-2j 



-U 



2/Q ,x(2n-l-2i)/2 



exp 



(36) 



by no more than a factor of o{u~ a ). Note that in the simplest case, when 
N = n = 1 and j = 0, this is precisely Lemma 2.4 of [2]. To prove the 
general claim we need here, observe that, by (35 ) and the assumption of the 
bounded support of the random vectors Uk, it is enough to prove that for 
every < j < L^^x^J an d 1 < m < n, 



(37) 



u 



**-l-2j]g . 



1 



-{n—m)/a 



-2/a 



m/2 



-U 



1 1 k 



exp 



o( u 



as u — > oo. To this end, note that the left hand side of (37) can be bounded 
by 



n-l-2j E 



2 L k 



2/a 



1/2 



2/a\( n - 2 i)/ 2 



exp 



(38) 



Take a small e > and write ( |38| ) as a sum of two terms, the first when the 

expectation is restricted the event {(J2k%2^ k 2 ^) 1 ^ > £u }-> an ^ the second 
arising when the expectation is restricted the complementary event. 

Noting that for any 6 > there is a finite c such that a~ 2S exp{— u 2 /a 2 } < 



cu 



-20 



for all u, a > 0, we see that the first term can be bounded by 



cu^E 



(£r-/") 1/2 i 



k=2 



V°° r" 



2/c 



1/2 



>eu 



Since P(Y,T=2 T k ^ > «) = o(n" Q/2 ) (see [8]), it follows that the first term 
in ( |38| ) is o(u~ a ) as u — > oo for every fixed e > 0. On the other hand, the 
second term is, clearly, bounded by 
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Since the sum T fc 2// " is a positive a/2-random variable, it has a density 
of asymptotic order x~( 1+a ^ 2 \ and straightforward estimates now show that 
the expectation above is of order u~^ n ~ 2 ^ +a ^ as u — > oo. Consequently, 



lim sup u a (the second term in (38)) < ce 

for some c > 0, and letting e —* proves (37). 

Returning now to ( 36 ) , which is what is left to study, we note that it has 
a simple structure, since it is immediate from the definition of W\ = W\{J) 
that it is a rank one matrix for all J G O n . Hence, 



det(r~ 2/a ^i 



I7 n/Q det (Wt) = 0, 



unless n = 1. Thus, taking n = 1 and j = (now the only possible value 
of j) in (36), recalling the independence of the ujk and T^, we need only 
consider N terms of the form 



E{|wi(i)|}E 



l/a 
■2/a 



1/2 



exp 



.N. 



The first expectation here is, by definition, /i,,, while the the second con- 
verges to the constant fioC a by the first part of Theorem 2.2 of [2]. Putting 
everything together proves (21). □ 



6. Concatenated-harmonisable fields 

For our final class of examples we shall introduce a class of random fields 
which, to the best of our knowledge, have not been studied earlier. We 
actually discovered them by looking for a class of examples which 'interpo- 
lated' between the sub-Gaussian ones, for which all the Cj appear in the 
asymptotic formula for the mean Euler characteristic of excursion sets, and 
the harmonisable ones, for which only Cq and C\ appear. However, having 
found them for this rather artificial purpose, we believe that they actually 
present an interesting class of stable fields that will provide useful models 
in applied settings. 

To define this new class of random fields, we take the representation 

/ C \ l/a °° N' 
m = (^\ ^ r,T 1/a E (Gu + G { $ sin(^)) . (39) 



Adler, Samorodmtsky and Taylor/Excursion sets of stable random fields 21 



where the are as in (19), the {Gy}, i = 1,2, I = 1, ...,N', are inde- 



pendent, standard Gaussian random variables, and the {ujkt} are indepen- 



dent with the distribution of the uj^ of (19). The parameter N' satisfies 



1 < N' < N. (When N' = 1 we recover the harmonisable fields of (19). 



We call such fields concatenated-harmonisable, the adjective "concatenated" 



coming from the innermost sum in ( 39 ) 



If the heuristics used before work again the dominant term in the expan- 



sion, 

1/a N' 



— — I r x 2^ \pkl cos{tu k e) + Gy smltuke. 



=i 

should determine the properties of the high level excursion sets. This random 
field is quite different from that of the simple random wave generated by 
the first term of the harmonisable processes, and so the arguments there, 
based on examples as in Figure [TJ are not going to carry over easily to 
the current situation. We did find an integral geometric argument which 



justified Theorem 6.1 below, but it was no longer simple and, since it was 
also non-rigorous, we shall not bring it here. 

In order to state the result, we need a little more notation. Changing 
slightly that of the proof of Theorem |5.1[ choose a facet J G O n , and, for 
each k > 1, define the n x n matrix Wk{J) with elements 

N' 

(Wfc(^))o-=E w M(*)wM(i), ( 4 o) 

for i,j € cr(J). Furthermore, define the (k independent) parameters 

A^EjldetO^J))! 1 / 2 }. (41) 
Theorem 6.1. Let f be a concatenated-harmonisable, SaS random field 



as in (39), defined on the N -rectangle T of (|2|). Assume that the control 
measure /j, has compact support, and a bounded density with respect to the 
Lebesgue measure. Then 

E{ip(Mf,T))} x u^^cJk^ + Y^KuY, \AKJ)],m 

V n=l JeO n ) 



where 

K 



2 



l-a/2 T /l±oA (jV')«/ 2 



irb a 
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and, for n = 1, . . . , N' , K n is given by (44) below 



If, furthermore, fi is rotationally invariant (so that f is isotropic) and M 
is a compact convex domain, then 



N 



E{cp(A u (f,M))} 



"VoC Q ^i^ n A„£ n (M), 



(43) 



n=0 



where A n is given by (41) for any n- dimensional facet J. 



Proof. We shall prove only (42), with the result (43) for isotropic pro- 



cesses on convex parameter sets following from the usual integral geometric 
argument via Hadwiger's theorem. 



The proof follows the lines of that of Theorem 5.1 and so we start by con- 



ditioning on the sequences {T^} and {uike}, to obtain a stationary Gaussian 
process on M. N with mean zero, variance 



- 2 = 7^V'E r , 



-2/a 



fe=l 



where 7 Q is as at (27), and with second order spectral moments 



k=i 



N' 

^2u ke (i)uJ ke (j). 
£=l 



We need to prove that the iterated expectation argument for computing 



IE {ip(A u (f, T))} is valid. This is done in Lemma 8.4 below. The arguments 
that worked in the harmonisable case also work here to show that the con- 



ditionally Gaussian process satisfies all the conditions of Theorem 3.1 



Thus, our task becomes one of evaluating (31 ) once again, albeit with the 



new definitions of the variables there. The term involving ^ changes from 
that in the previous proof only insofar as there is now an additional factor 



of N' in the definition of a in (28) and this gives the first term in (42) (cf. 



Lemma 4.2 



As far as the other terms are concerned, the argument is identical to that 



in the proof of Theorem 5.1 as far as (36), and so what remains to compute 
is 



^ N fy(2n-l-2j)/2 E 




exp 



2/a 
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for 1 < n < N' and < j < L^r^J • The restrictions on n follow from the 
fact that the matrices Wk{J) have rank not exceeding N' and, since the Ukt 
have a distribution with a non- vanishing absolutely continuous component, 
this rank will equal min(n, N'), with positive probability. 

Since W\ is independent of the r&, we can rewrite the above as 



(iV / )- (2n - 1 - 2j)/2 E{|det(W 1 (J))| 1 / 2 } 
xE < 



p-«/ Q n n-l-2j 



-U 



2/a N (2n-l-2j)/2 



exp 



2/a 



The first expectation here is, by definition, A(J). As for the second, its 



asymptotics follow from Lemma 6.2 below, to give 

u n - l - 2 i x ^C a C nj u-^ n - l -W = fi C a C nj u- a , 

where 

Cnj 4 aC Q/2 </V2 (a+n - 2j - 3)/2 r ( a+ "- 1 - 2j ) 7 ™-i-2i (iV ')(«-n)/2. 



If we now substitute this back into (31 ), collect all the constants appearing 



in the Hermite polynomials, and define 

I — I 



-iv a 



nj 



( 27r )(n+i)/2 ^ j!( n -l-2j)!2i 



(44) 



a few lines of algebra yield (42), and we are done. 



□ 



Lemma 6.2. Maintaining the above notation, set 



X = Y / TJ 2/a , n>0, p>n/2, 

3=1 



and take 7 > 0. Then 

E ^ T -n/a x -f3 e - u y 2 ^X^ 



(45) 



a2{a+ 2^n~2)/2 Ca/2< /2 ia+ 2p-n T{(3 + 1 (q _ ^^-(a+2/3-n) _ 
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Proof. The proof will proceed by establishing asymptotic upper and lower 
bounds for the expectation, which we shall denote by l u . 

For the upper bound, note that since T 1 2 ^ a < X it is immediate that 

/ u <E{x-^-"/ 2 )e-" 2 /27 2 X}_ (46) 



However, the asymptotics of the right hand side are covered by Lemma 4.2 



and so we have that the right hand side of (45) provides an upper bound 
for the asymptotics of I u . 

For the lower bound, fix e > and again exploit the fact that T 1 2 ^ a < X 
to see that 

> (1 + e)-% jr-/"^-^-/^ wrr 2/Q ^ ^ a/B<a , _ a/a } 
where Y is a copy of X independent of T\. Now fix M > and note 

iu > (i + e)^F{x < m}e {r-^Crr 270 )-^-" 2 /^ 2 ^ 270 !^,/^^} 

= (1 + e)-^{X < M}E |(r- 2/a )-^-/ 2 ) e -" 2 / 2 ^ 2r " /a l r „ 2/Q>M/ J . 



4.2 



with r x 2 ^ a 



The expectation here is similar to that in (46) and Lemma 
replacing X. However, these two random variables have precisely the same 
tail behavior, and a check of the proof of Lemma |4.2| (cf. Lemma 2.2 of 
m) shows that this is all that entered into the asymptotic behavior of the 
expectation. Hence 

Iu > (i + £ )-P F {x < M}E jx-^+^e-" 2 / 2 ^} . 
Sending e — > and M — > cxd completes the proof. □ 



7. On the mean Lipschitz-Killing curvatures of excursion sets 



Throughout this paper, we have concentrated on the expected Euler char- 
acteristics of excursion sets. However, at least in the isotropic cases, these 
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results are immediately extendable to expected Lipschitz-Killing curvatures, 
via a result known as Crofton's formula. 

To state Crofton's formula we start with the affme Grassmanian Graff (TV, k) 
of all fc-dimensional flats in M, N ; viz. of all /c-dimensional linear subspaces 
of M. N not necessarily passing through the origin. On Graff (JV, k) there is 
a natural Haar measure, known as kinematic measure, which we denote by 
X^. Its precise definition and normalisation will not be important to us. 
Crofton's formula states that 



CJMnv) dx%_ k (v) 



GiaS(N,N-k) 



k + j 

3 



C k+j {M). 



(47) 



Now let / be one of the stable random fields of this paper, and M a 



compact, convex set in R . If / is isotropic, then, recalling that the Euler 
functional (p is also Co, all of our results can be written in the form 



N 



E{C (A u (f,M))} = J2C k (u)C k (M), 

k=0 



(48) 



for some functions C k {u). The C k {u) depend on the parameters of / and, 
quite often, are identically zero. What is important, however, is that they 
are dependent neither on M nor on N = dim(Af). 

Since if M is compact and convex so is M n V, for any V £ Graff (AT, k). 



Thus, exploiting Crofton's formula twice, and (48) once, we have 



E{£j(4(/,M))} 



E<! / C (A u (f,M)nV) d\%JV) 

'GraS(N,N-j) 



E {Co (A u (f, M) n V)} d\%AVl49) 



GraE(N,N-j) 



N 



Graff(/V,7V-j) ^ 



N 

E C k {u) 
k=o ( 

N—j 

E 

fc=0 



j + k 
k 



Gra,S{N,N-j) 

C k {u)C j+k {M) 



C k {MDV) d\%_j(V) 



the change on the range of summation coming from the fact that Ci{M) = 
for all I > N. 



Adler, Samorodnitsky and Taylor/ Excursion sets of stable random fields 26 



Using this, all the formulae in this paper for the asymptotics of mean 
Euler characteristics, in the isotropic cases, can be extended to asymptotics 
for mean Lipschitz-Killing curvatures. 

Of course, something has to be said about justifying the change of order of 



integration and expectation at (49), which we shall do in Lemma 8.5 below. 



8. Appendix: On the finiteness of expectations 

A serious technical point which we avoided throughout the paper was justi- 
fying the conditional expectation arguments 

E{<p(A u (f,T))} = E{E{<p(A u (f,T)) (50) 

where / was stable, and T was the information we assumed to make / 
conditionally Gaussian. 

The need for justification lies in the fact that Euler characteristics need 



not be positive. (Were they positive, (50 ) would always hold, with both sides 



being finite or infinite together.) We shall show that 

E{\<p(A u (f,T))\}<oo, (51) 



which is sufficient for (|50j), restricting ourselves to rectangular parameter 
spaces only. The same arguments can also be applied for general compact 
domains, the only difference being a heavier investment in notation. 



Before starting the proof of (51 ), we note that all the random fields that 
we considered in this paper are suitably regular Morse functions, in the 
terminology of Chapter 6 of [3 . This follows from a Fubini argument and 
the fact that, with probability one the conditionally Gaussian random fields 
on which we based all our calculations are of this kind. This, in turn, follows 



from the fact that they all satisfy the conditions of Theorem 3.1 which, by 



Corollary 11.3.2 of [3 j5 implies that they are suitably regular Morse functions. 



Turning now to the proof of (51 ), let J, as usual, be a n-dimensional facet 
of a iV-dimensional rectangle T . Then, for a Morse function, the Euler char- 
acteristic of the excursion set A u (f,T) can be represented as an alternating 
sum, over all J, of the numbers of critical points of various indices of f\j 
over A u (f\j, J) (cf. Section 9.4 of [3]). This sum is trivially no greater, in 
absolute value, than the total number, over all J, of the number of critical 



points of f\j over J. Thus, in order to establish (51), it will suffice to show 
that for every facet J 

E{E{iV(/,J)|^}} <oo, 
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where N(f, J) is the number of critical points of / when restricted to J. 

In order to compute the inner expectation, recall / is conditionally a 
stationary Gaussian random field in all the cases of interest to us. Call this 
process /. Then it follows immediately from Theorem 11.2.1 of |3] and the 
independence of first and second order partial derivatives of Gaussian fields 
(cf. [3] Section 11.7) that 



E 



EjldetV 2 ^)} 



detS 



(52) 



f., 



where V 2 /|j is the n x n matrix of second order derivatives of /|j and S- 

the n x n covariance matrix of its first order derivatives. 

With these preliminaries behind us, we can begin establishing the two 



main results of this appendix, Lemmas 8.1 and |8.4| 

Lemma 8.1. Let f be a sub-Gaussian process satisfying the conditions of 



Theorem 4-1 Then, for every ii£l, K{\tp(A u (f , T))\} is finite. 



Proof. Recall that in the sub-Gaussian case / is no more that X l / 2 g, 
where g is Gaussian and X is positive a/2-stable. Thus, given X = x, the 
conditioned process / of (52 ) is no more that x l l 2 g. It therefore follows that 



for a facet J of dimension n, 
E{N(f,J)} = 



e{e {iV(/,J)|x}} 
E{|detV 2 ^|} 



E 



E 



detS 



fj 



X 



X n / 2 E{|detV 2 5 |j|} 



X"/ 2 detS 



91 J 



X 



E{|detV 2 5 |j|} 



detS 



g\j 



The last expression is purely Gaussian, and under the assumed conditions 
of non-degeneracy, also clearly finite, so we are done. □ 



Before turning to concatenated-harmonisable processes we require a tech- 
nical lemma. 
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Lemma 8.2. Let A and B be two n x n matrices, and suppose that A is of 
rank m < n. Then each term in the standard Laplace expansion ofdet(A+B) 
involves at most m elements of A. 

Proof. We shall prove the result for m = 1. The result for general m then 
follows by writing a matrix of rank m as the sum of m matrices of rank one. 

Since A has rank one, there is a vector a and numbers 6\, . . . ,6 n such that 
the j-th row of A is 6ja, j = 1, . . . , n. Let bj be the j'-th row of B. Writing 
the determinant of A + B as a function of its rows we have 

det(A + B) = det(9 1 a + b 1 ,9 2 a + b 2 ,...,9 n a + b n ) 
= 6>idet(a, 9 2 a + b 2 , . . ■ , 6 n a + b n ) 

+det(6i, 9 2 a + b 2 , . . . , 9 n a + b n ). 

Note that 

9±det(a, 9 2 a + b 2 , . . . ,9 n a + b n ) 

= 9\9 2 det(a,a, . . . ,9 n a + b n ) + 6>idet(a, b 2 , . . . , 9 n a + b n ) 
= 6 , idet(a, b 2 , . . . , 9 n a + b n ) 
= ... = 6»idet(a,6 2 , b 3 ...,b n ). 

Similarly, 



det(6i, 9 2 a + b 2 ,...,9 n a + b n ) = ^ # fc det(&i, . . . , a, b k+1 , ...,b n ), 

k=2 

and continuing this process leads to 

n 

det(A + B) = detB + ^ 9 k det(bi, . . . , b k - 1: a, b k+ i, . . . ,b n ), 

k=l 

from which the result follows. 



□ 



An immediate consequence of Lemma 8.2 



is 



Corollary 8.3. Let Ai, . . . , A n be n x n matrices of rank m < n. Then, for 
any real r±, . . . , r n , 



r l . . . r n 



3=1 
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We now have what we need to prove our last lemma. 
Lemma 8.4. Let f be a harmonisable, or concatenated-harmonisable, SaS 



random field satisfying the conditions of Theorem 5.1, or Theorem 6. 1 
spectively. Then K{\ip(A u (f,T))\} < oo, for every u G R. 



re- 



Proof. We shall prove the result in the notation of the concatenated- 
harmonisable case, since taking N' = N gives the harmonisable case. In 
view of the argument leading to ( 52 ) , and the st ruct ure of the conditional 
process / in this case (cf. the proof of Theorem 5.1 for this and following 
notation) we need to show that 



E 



E 



{|det£r=i^ 1/Q GW fe |} 



detET =1 r~ k 2/a w k 



1/2 



< OO, 



(53) 



where the inner expectation is taken only with respect to the Gaussian 
random variables. Recall also that each W k , given by (40), can be written 
in the form 

N' 

W k = Y,W k ,l, (54) 
1=1 

where (W k ,i, k > 1, I = 1, . . . , N') are i.i.d., rank 1, n x n, random matrices 
of the form ( 33 ) . 



We concentrate on the denominator in (53) first. Set 

,1 H 



N' 



and recall that the determinant of the sum of non-negative definite matrices 
is at least as large as the determinant of each of the terms, so that 

/ oo \ 
-2/at, 



det £ r r^ 



\k=l 



>det (j:r k 2/a W k + T~l{ a (w d+ 

\k=l 



11 + ... Wd+l, n -dN' 



-2/a r -2/a 
: 1 d+1 



\k=l 



N ' ( " 

| det W k + (Wd+i,! + . . . W d+l 

,n-dN> 



\k=l 



the last line following from Corollary 8.3 Set 

/ d 

D = det [J2 W k + (W d+1>1 + ... W d+1 

\k=l 



n-dN' , 
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so that, to prove (53), we need only check that 



E 



E{|det£^ =1 r fe 1/Q G fc Ty fc |} 



D (TV 1 r -2/a r -2/a 
^(llfc=l i fc L d+1 



N' 



< OO. 



(55) 



For later reference, note that D is a polynomial function of n i.i.d. ran- 
dom vectors, each of length n. Under the conditions of the theorem, the n 2 
random variables defining D have a bounded joint density with compact 

2 

support in W 1 . We claim that it follows from this that there exists a finite 
C such that, for e small enough, 



P{£> <e}< Ce. 
In turn, from this it trivially follows that 

E{ZT 1/2 } < oo. 



(56) 



(57) 



To prove (56), first write X\, . . . , X n for the n random vectors defining the 



determinant D, and let A C 1" denote the support of each one. Fix e > 
and define 



A, 



(xi, . . . ,x n ) e A n : min \xi(j)\ < £ L 

l<t,j<n ) 



Then 



B £ = {x = (xi, . . . , x n ) G A n : D(xi, . . . , x n ) < e] . 



P{D <e}<v (A E ) + v(B e n (A n \ A £ )) 



(58) 



where v is the measure on A n generated by the random variables Xi, . . . , X n . 
Both of these terms are easily seen to be bounded by a constant multiple of 
e, of e is small enough. The first is small, since A e can be covered by a finite 
number of n-dimensional rectangles, one of whose sides has length no more 
than 2e and the remaining sides having length no more than the diameter of 
A. Since v has bounded density, this gives that \i (A e ) < Ce for some finite 
C. 



As far as the second term in (58) is concerned, we note that, because of 
the smoothness of the mapping D : W 1 — > M in the region A n \ A £ , the set 
B e D {A n \ A £ ) is a C°°, locally convex, stratified manifold, and so it follows 
from the generalized tube formula of |3 (Theorem 10.9.5) that this set also 
has i/-measure bounded by Ce. 
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Returning now to ( 55 ) , note that by expanding the determinant and ap- 
plying Lemma 8.2 we have that (55) is equivalent to 



E ( 



E {|E£=i • • • EZ=i T ~kl' a ■ ■ ■ <' a G^ ■■■Gk, B kl ... kn l F \} 



D 



(n 



r -2/a p -2/a 
fe=l 1 k 1 d+1 



where the B^..,^ are uniformly bounded random variables, independent of 
Tj and the Gaussian Gj , and F is the event that at most N' of k\ , . . . , k n 
are equal. 

Fixing the Tj, applying the boundedness of the B^...^, taking an expec- 
tation over the Gaussian Gj and using their symmetry, bounds the above 
by 



CE 




-2/a 



-2/a 



1/2' 



ki=l 



(ntirr'i 



2/a r -2/o 
k 1 d+1 



(59) 



where C is a finite constant which may change from line to line. Now use 
the fact that D is independent of all the other random variables in the 
above expectation, along with (57), to remove the factor of D^ 1 from the 
expectation, with an appropriate change of the constant C. 



We may, and shall, assume that, in the fe-fold sum in (59), k\ < k-z < 
• • • < k n . Thus what remains to show is that the remaining expectation, 
which we rewrite as 




■ E 

fc n _l=fe„_2 + l 



p-2/a p-2/o 



(ni 



i l k 



,-2/q 



N' 



(r^°) 

oo 

E 



iV'-l 



1/2 



T -,2/a-p-2/a 1 
1 d+l 1 fe n l F 



L+l 



is finite. Consider the summation over k n , keeping ki, . . . , k n _\ fixed. Since 
each term in the (n-fold) sum is bounded from above by 1, we can begin the 
summation over k n at an arbitrary large value K. Then, the expectation in 
the last sum is bounded by 



1/2' 



E 



rti E ?l 



(60) 



-K 
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Since Gamma random variables have all moments finite, by Holder's inequal- 
ity it is enough to prove that for e < 1, 



K k n =K 



> < OO. 



(61) 



Furthermore, since it is easy to check that 

E{r fc 2/Q } <ck~ 2 ' a 



for k large enough, the finiteness of (60) now follows from the fact that the 

(l-HO/2 



expectation in (61) is bounded by 

oo 

E 



,k„=K 
and that a < 2. 

Next, we consider the double sums obtained by fixing ki, . . . , k n -2 and 
taking the sum only over /c n _i and k n . As in the case we have just considered, 
we can start the sum over k n ^\ at a (large) K of our choice, and then the 
expectation involving the double summation, with an additional factor of 
^ 2 J+i brought into play, is bounded above by 



oo oo 



1/2' 



MK+i E ECS 

k n _i=K k„=K 



The finiteness of this expression follows as above, from the easily checkable 
fact that 

E{r^ /a r^ /a }<ck- 2/a k 2 2/a 

for all k\ < &2 large enough. Iterating this argument, we arrive at the finite- 
ness of (59), once we prove that, for all k\ < &2 < • • • < k n large enough, 



-2/a 



This, however, follows from Holder's inequality and the fact that for all k 
large enough, 



E{r k 2n/a } <ck- 



-2n/a 
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This completes the argument. 



□ 



Lemma 8.5. The exchange of expectation and integration in (49) is justi- 
fied. 



Proof. In the notation of (49), it will suffice to show that 



Giaff (N,N-j) 



E{|£ (A u (f,M)nV\)} d\%^{V) <oo. 



(62) 



Note first that ( 49 ) was only being applied under the conditions of one of 



the main theorems of the paper, so we know that the inner expectation is 
always finite. We also know, from the discussions of this appendix that, for 
any convex A, 

|£ (A u (f,A))\<^(A), 

where we define ip{A) to be the number of critical points of / in A and 
on its various boundaries. All of these also have finite expectations. The 
functional IE{V ; } is clearly additive (since ip itself is) in the sense of (§, as 
well as being invariant under rigid motions (by isotropy) and continuous 
in the Hausdorff metric. Thus, by Hadwiger's theorem (cf. (|7|) there exist 
constants Cj, dependent of u and the structure of /, such that 



N 



EM 



c k Ck(A). 

k=0 



Substituting into ( 62 ) we therefore have 



L 



E{|A) (A u (f,M)nV)\} d\%_AV) 



Gra,S(N,N-j) 



< I EWMny)} d\%_AV) 

lGiaS{N,N-j) 

N 

Y / c k £ k (MnV)d\%JV) 

lGiaff(N,N-j) ^ 
N 

E C W £ k (MnV)d\%_ j (V) 

^ JGr a S(N,N-j) 
' \k + 7I 

' ' C k+j (M), 



3=0 



k 
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the last line following from Crofton's formula (47). Since M is compact and 
convex, the Cj(M) are all finite, and so we are done. □ 



References 

[1] Adler, R. and Samorodnitsky, G. (1997). Level crossings of abso- 
lutely continuous stationary symmetric a-stable processes. Ann. Appl. 
Probab. 7, 460-493. 

[2] Adler, R., Samorodnitsky, G. and Gadrich, T. (1993). The ex- 
pected number of level crossings for stationary, harmonisable, symmetric, 
stable processes. Ann. Appl. Probab. 3, 553-575. 

[3] Adler, R. and Taylor, J. (2007). Random Fields and Geometry. 
Springer Monographs in Mathematics. Springer, New York. 

[4] Adler, R., Taylor, J. and Worsley, K. Random Fields and Ge- 
ometry: Applications. In preparation. 

[5] Breiman, L. (1965). On some limit theorems similar to the arc-sine 
law. Theory of Probab. Appl. 10, 323-331. 

[6] Hadwiger, H. (1957). Vorlesiingen Uber Inhalt, Oberflache und 
Isoperimetrie. Springer-Verlag, Berlin. 

[7] Klain, D. and Rota, G.-C. (1997). Introduction to Geometric Prob- 
ability. Cambridge University Press, Cambridge. 

[8] Samorodnitsky, G. and Taqqu, M. (1994). Stable Non-Gaussian 
Random Processes. Chapman and Hall, New York. 

[9] Worsley, K. (1997). The geometry of random images. Chance 9, 
27-40. 

Robert J. Adler 
Industrial Engineering and Management 
Technion, Haifa, Israel 32000 
e-mail: robert@ieadler . technion .ac.il 
url: ie . technion . ac . il/Adler . phtml 

Gennady Samorodnitsky 
School of Operations Research 
and Information Engineering 
Cornell University, Ithaca, NY, 14853. 
e-mail: gennadyOorie . Cornell . edu 
url: www . or ie . Cornell . edu/ ~gennady/ 

Jonathan E. Taylor 
Department of Statistics 
Stanford University 
Stanford, CA, 94305-4065. 



Adler, Samorodnitsky and Taylor/Excursion sets of stable random fields 



35 



e-mail: taylor@dms . umontreal . ca 



url: www-stat . stanford.edu/~jtaylo/ 



